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1. INTRODUZIONE 


Nel seguito verranno esposti risultati di esistenza ed uni- 


cità relativi al sistema di equazioni ordinarie 


(1) AX + Bx = f 


Tini indica l'insieme delle matrici su C con h righe e 


ove A, BeC 
k colonne: è usuale - in gran parte della letteratura sull'argomento - 
svolgere la teoria di (1) in riferimento alle matrici, sebbene questa 

sia essenzialmente relativa agli operatori canonicamente identificati 

ad esse), f (il termine noto) è una funzione R + C" e 1a incognita 


(1) 


t + x(t) viene ricercata in € ec è chiaro che (1) è interes- 

sante se A non è invertibile. 

Nonostante sia relativamente elementare, la (1)merita già 

un certo interesse per vari motivi: 

- sebbene la teoria sia stata elaborata in anni recenti ( citiamo in 
particolare S.L. Campbell) essa ha già raggiunto risultati conclusi, 
che sono estensioni, per molti versi completamente soddisfacenti, dei 
risultati e delle tecniche classici (caso di A invertibile); 

- al tempo stesso, già nello studio di (1) emergono gli strumenti (in 
particolare, un tipo generalizzato di inversa per A) che sono i me- 
glio idonei a trattare casi più complessi; 

- la possibilità di studiare, sulla base di una esauriente conoscenza 
di (1), casi più complessi di sistemi degeneri, con le usuali tecni- 
che di estensione: ad esempio, il caso di coefficienti dipendenti dal 
tempo ([3]), il caso della dimensione infinita ([1], [4]) - per cui tut 


tavia la teoria è lontana dall'essere approdata a risultati conclusi- 
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vi-,il caso di equazioni semilineari, il caso dei sistemi d'ordine su 
periore ([3]); 

- il numero notevole di applicazioni di (1) (alcune sono riportate nel 
testo e nella bibliografia di [3]), che in molti casi hanno precedu- 
to la teoria completa, indice questo sicuro di interesse applicativo. 

I] testo seguente è diviso in due paragrafi: nel primo vengo- 
no rapidamente esposte le principali proprietà dell'inversa di Drazin 

di una matrice, che, nell'altro paragrafo, è il principale strumento 


per studiare la (1). 


2. L'INVERSIONE SECONDO DRAZIN 


2.1. Per motivare l'introduzione di un nuovo tipo di inversa 
nel trattare (1), ricordiamo che nel caso 'classico' in cui esiste n°! 
il problema di integrare (1) è essenzialmente algebrico o, più esatta- 


mente, di teoria spettrale (in senso lato) delle matrici; si tratta di 


calcolare expla 80); poi, moltiplicata la (1) per exp(a '80)a”!, ci 
si riduce a 
i (exp(a '8t)x) = exp(a° ‘Bia !r 
risolvibile per integrazione: 
-1 la -1 x 
(2) x(t) = exp(-A_ Bt)C :f exp(-A B(t-s))A f(s)ds 
lo) 


n è : 
ove C = x(o)e C è arbitraria. 
Questa classica formula suggerisce due considerazioni, utili 


all'estensione di questi risultati: 
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a) nella (2) compare in maniera essenziale Pao accanto all'esponenzia- 
le, determinabile concretamente con l'aiuto di proprietà spettrali 
dele matrici considerate: come si vedrà, l'inversa di Drazin, a dif 
ferenza di altri tipi di inversa generalizzata, mantiene proprietà 
siffatte il più vicine possibile a quelle dell'inversione ordinaria; 

b) la somma, in (2), riflette il ben noto principio di sovrapposizione 
che vale anche nel caso degenere: anche qui verrà studiato anzitut- 


to il sistema omogeneo associato a (1), 
(3) AX + Bx = 0. 


2.2, Come già accennato, i vari processi di inversione gene- 
ralizzata che consentono di trattare il caso di sistemi degeneri di 
equazioni di 1° grado, non mantengono in genere le proprietà spettra- 
li dell'inversione usuale, e quelle correlate: ad esempio, se A è in- 


vertibile, gli spettri o(A) e #4 sono legati da 
-1 
(o(A)) = o(A 


_ - : (9 4 
relazione che non vale în genere per autovalori # 0 se ad A si sosti- 


tuisce l'inversa di Moore-Penrose aÎ di A (se 


1 -] po 1/2 0 4 
A = | è A | e si ha subito o(A) = {0,1}, 


0 0 =1/2 0 
o(A°) = {0,1/2}); questo tipo di proprietà si estende meglio al caso 
singolare con il tipo di inversione proposto, in un ambiente algebrico 
più astratto (anelli) da Drazin nel 1958 ([61): per introdurre la de- 


finizione è però opportuno premetterle i] 


Teorema 1. ([2, p. 1211). SeAe pre vi è un intero k > 0 


(detto L'indice di A, ind (A)) che è il minimo intero = 0 tale che, 
equivalentemente: 


i) N (AS) «UO, oppure 
ii) R(AS) = R(A°4%); 
k 


per questo k vale La decomposizione cla RA“) MA). 


(R(A), N(A) indicano l'immagine ed il nucleo di A; © indica come sem- 
pre somma diretta). 

11 numero ind (A) è una misura della non .invertibilità di A, perlome- 
no nel senso che ind(A) = 0 @ A è invertibile (si conviene che per B 


matrice nxn, sia 8° = identità): sebbene enunciato per ma- 


trici, il teorema concerne evidentemente l'operatore canonicamente as- 
sociato ad A, cui si riferisce anche la nozione di indice; con questa 
osservazione è chiaro come la seguente definizione sia anch'essa in- 


trinseca. 


Definizione 1. Stano Ae CU", k = ind(A). SÉ dice inversa di 
Drazin di A una matrice Pa tate che 


i) Pad; 


In effetti, se si indebolisce iii in iîi' : es<ste un pe N 


Fa AP, si ottiene una definizione equivalente a questa: 


per cui ap! A 
anzi, l'originaria definizione di [6] consta proprio di i-ii-iii'. 


La correttezza della definizione si basa sul seguente 
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Teorema 2. Ogni matrice A € C" ammette un'unica) inversa di 
Drazin n, 

Per la sua natura costruttiva, riportiamo le grandi linee 
della dimostrazione (cfr. [2]), marcandole con. 

mCon riferimento al teorema 1, è facile vedere che R(A“), vas) 

sono sottospazi invarianti per A, la quale sul primo è invertibile men- 
tre è nilpotente (d'indice k) sul secondo. Per la definizione 1, questi 
sottospazi sono invarianti per ogni eventuale matrice Mu del tipo detto: 
bisogna che p su R(AS) inverta A e, su N(AS), che sia a = 0: questi 
requisiti individuano univocamente A.m 

Osservazioni.l. Se k = 0, cioè A è invertibile, è n = 1, 
se k > 0 la decomposizione in somma diretta del teorema 1 fornisce per 


A la rappresentazione 


(4) A=P P 


xn _, ui dia 
con Pe c" invertibile, C quadrata invertibile, N quadrata nilpotente 


d'indice k, sicché in corrispondenza di (4) 1a costruzione del teorema 2 


che fornisce una rappresentazione a blocchi per a 


2. Per come è definita, AD non è adatta alla risoluzione di 


I- 8. 
Ax = b, cioè, in generale, ab non è soluzione di questo sistema, an- 
che se è b € R(A): ciò avverrebbe se fosse aa = A (se b = Ac, allora 
A(APb) = marc = Ac = b), ma si ha subito che ana = Asee solo se 
ind As 1. 
m Se vale ma = A ed ind(A) = 1, con riferimento a (4)-(5) 


si ha 
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da cui N = 0, cioè ind A = 1; viceversa, trascurando il caso banale 
ind(A) 


della definizione 1. 


0, se ind(A) = 1, l'equazione da provare coincide con 1a iii 


Quindi lu è un'inversa generalizzata algebrica nel senso di 
Nashed e Votruba ("A Unified Operator Theory of Generalized Inverses", 
in [8]) se e solo se ind(A) s 1: in tal caso essa viene in genere chia 


mata inversa di gruppo (cfr. [2]). 

2.3. Come si è detto, l'inversione secondo Drazin permette di 
mantenere per Pa diverse proprietà - specie quelle connesse alle spet- 
trali -, tipiche di "a nel seguito menzioniamo alcune di esse, riman- 


dando per la dimostrazione di queste e di altre, a [2]. 


Teorema 3. Se A € ig e si pone, per VE C 


A sex 0 
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nisulta XE o(A) © x e o(A°). 

Se X € o(A), X # 0, x è un autovettone generalizzato di grado 
p per A relativo a X (cioè (A-x1)Px = 0, x # 0), 4e e solo se x è un 
autovettone generalizzato di grado p per Pa nelativo a xÎ = i, 

La dimostrazione è in sostanza una verifica, sulla base di 
(4)-(5): per un inquadramento più generale di questo risultato, si ve- 
da Ben-Israel e Greville, "Some Topics in Generalized Inverses of Ma- 
trices" in [8]. 

Un'altra proprietà di cui a, a differenza di altre inverse 
generalizzate, gode sempre, e che ha rilevanza anche per il calcolo ef- 


fettivo di n, è la seguente: 


Teorema 4. Se A € pi c'è sempre un polinomio p a coeffi- 
cienti complessi tale che n° = p(A). 

Come applicazione del teorema 4, consideriamo una proprietà 
dell'inversione che non si conserva in generale per alcun tipo di in- 


tri, essa 


versione generalizzata, la regola del prodotto (18)! = BO 
non vale nemmeno per l'inversa secondo Drazin, anche se in questo caso 
la formula può essere stabilita sotto la sola ipotesi AB = BA, a diffe- 


renza di altri tipi di inversione: 


Teorema 5. Se A, BE riti commutano , 44 ha 


(a8)? = n°g? Li pe, 


tm Con la rappresentazione del teorema 4, è facile verificare che 
D D 
A,A ,B,B commutano; si può allora vedere facilmente che ndg? = 8PpD ve- 
rifica i-ii-iii' nella definizione 1 per (18). 


a RES 4 4 è D 
Menzioniamo infine una formula di approssimazione per A . 


Teorema 6. Se A € ca allona se L= ind(A), 


PETE alt 1 UP 


2>0 
Inoltre si può provare che 0 < ind(A) = k è anche l'ordine della sin- 


golarità polare di x + R(A,A) nell'origine, nel qual caso si ha 


lin! R(z,A) = canP_ 1a! 
Z>0 


(cfr. anche [7]). 


2.4. Sono state proposte diverse estensioni - sinora non comple 
tamente soddisfacenti - dell'inversione secondo Drazin in spazi di Ba- 
nach qualunque. Nel caso in cui T sia un operatore a indice (nel senso 
che a.(T) = S(T) < +: cfr. [5]) con polo di ordine k = ind(T) per 


z + R(z,T) nell'origine, si può porre 


t. tea! i x ROTA 
|a|=e 


e sviluppare una teoria assai simile a quella in dimensione finita: tut- 
tavia questo impone notevoli restrizioni alle applicazioni. Diverse al- 
tre varianti sono state proposte (vedi in particolare Campbell, "The 


Drazin Inverse of an Operator" in [4]). 


3. L'EQUAZIONE AX + Bx = f » A SINGOLARE 


3.1. Cominciamo col considerare alcuni semplici esempi di (1). 
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010 001 
A) Siano A = [o 0 .| , B= [o 0 o| » f= 0: (1) diviene 
000 00 


0; 


x 
U1] 


0 
si hanno soluzioni solo con dato iniziale x(0) = fa) ed in tali casi 


s } A B 
la soluzione non è mai unica. 


100 
B) Siano A = [o 0 | » B= I: allora il sistema (1) è 
000 


x +% ah, 
Sta 
e ia 


baia 116 ANPIROA a 
ciò richiede f.e C°‘: se ciò è vero, Xo = fo 3 i * % + f, dà 
soluzione unica se si impone una condizione iniziale del tipo 


o . 
x(0) = £(0) - î,0| » 0€ C arbitrario. 
£,(0) 


Dunque. nel caso che A sia singolare: 


1) il dato iniziale è soggetto a vincoli; 
2) il termine noto f nella non. omogenea deve essere, in generale, abba- 


Stanza regolare, se si vuole che il problema di Cauchy 


AX + Bx = f 
x(0) = x € È 


(6) 


n 


abbia soluzione; 
3) anche nelle ipotesi 2 la soluzione può non essere unica. 

Nel seguito tratteremo il caso del problema (6) per la (1) 
considerando solo soluzioni ‘classiche’ x € ct!) (R, c"): del resto in 
dimensione finita le stesse ipotesi che garantiscono esistenza ed uni- 
cità implicano la possibilità di prolungare soluzioni locali (in dimen- 
sione infinita le cose stanno ben altrimenti: per tutta la questione 
cfr. Campbell cit., in [4]). Non stupisce dunque i1 prevalere dell'aspet 


to algebrico nella trattazione di (6), come nel caso di A invertibile. 
3.2. Introduciamo la seguente terminologia: 
Definizione 2. x si dice un vettore (0 condizione) iniziale 
compatibile con (1) se il problema (6) ha (almeno una) soluzione; La 
(1) si dice trattabile se (6) ha soluzione unica per ogni condizione int 
ziale compatibile. 
Nel seguito considereremo soprattutto la trattabilità di (1): 


seguendo un procedimento euristico classico, si può applicare în (6) 


una trasformazione (formale) di Laplace: 
L(A&X) + L(Bx) = L(f) 

da cui 
(SA + B)L(x) = L(f) + Ax 


che suggerisce una stretta connessione tra unicità ed esistenza del ri- 


solvente generalizzato s + (SA + 8); in effetti vale il 


Teorema 7. ([3, p. 34]). Se A, BE o" L'equazione omogenea 


(3) è trattabile se e solo se: 

(7) 3X€ C tale che XA + B è invertibile . 

Poiché questo teorema è la chiave della trattazione di (1) in [3], dia- 
mo le grandi linee della dimostrazione . 

mValga (7): definiamo le matrici À = (2A + 8)" lA,8=(24+8)'8; 
I 


poiché AA + B = 4 À e B commutano e (3) equivale a Ax + Bx = 0; a me- 


no di similitudini si può scrivere 


Cc 0 
A = | ] con C invertibile, N nilpotente, da cui 
ON 


_ I- AC, 0 Bj 0 x 
s-[ La | | ; posto x= [1] 
0 I-AN 0 B 


9 2 2 


(3) si spezza nel sistema 


" 
(| 


Ck, + (I-1C)x, 


" 
(©) 
we 


Nk, + (I-xN)x, 


poiché la prima equazione è chiaramente trattabile, basterà provare che 


tale è la seconda; a questo scopo proviamo che ha la sola soluzione x.= 0. 


Fa 
In effetti, se k = ind(N), applicando NÉ! segue ((I - N) è invertibi- 
le) nix, = 0; quindi niet, = 0 e, applicando pra. nie, = 0 e co- 


sì via. 


Viceversa, (3) sia trattabile e, per assurdo, 


veca ve e, we 
tale che (XA + B) x, * 0. 


Si proverà che in tale ipotesi si può costruire una soluzione non nul- 
la del problema (3) con condizione x(0) = 0, ciò che contraddice la trat 
tabilità di (3). 

Una tale soluzione si costruisce così: se 
Arese) (RAS ++ SRO) sono tali che i corrispondenti v, sono linear- 


mente dipendenti, si scelgano scalari non tutti nulli dj +9,, per cui 


r 
* a. v. = 0; allora z(t) = = a, exp(i.t)v. è la 


soluzione cercata. Ml 
Dal teorema 7 segue una caratterizzazione completa della con- 


sistenza e della trattabilità per l'omogenea (3). 


Teorema 8. ([3]). Se (3) è trattabile, cioè vale (7), allona 
i) un x € c" è consistente per (3) se e solo se 

LE RATA); 

Co) 
ji) L'‘ntegnale generale di (3) è dato da 


(8) siti = pil Pia, dec. 


Nell'enunciato si sono seguite le notazioni introdotte nella prova del 
“De Da 4 
teorema 7: si può peraltro dimostrare che A A, A B non dipendono da A. 


Si noti ancora l'analogia formale tra (8) e la (2) nel caso omogeneo. 


I- 15, 


3.3. Per l'equazione non omogenea (1), oltre alla {7), bisogna, 


come si è visto, fare alcune ipotesi sulla regolarità di f. Si ha 


Teorema 9. ([3]). Sé consideri L'equazione (1), nell'ipotesi 
(7); di più, posto k = ind(A), sta f € cl (rc). Allora La (1) ha sem- 
pre soluzione,una:essendo data da 


È 
x(t) = cap(-A%) [ exp(Ags) Rf (s)ds + 
A i 


POR n'a lt) rt). 
i=0 


Posto W = r(0), 4% ha: 
& -D- 
i) x € c" è consistente per (1) © x € W + R(AÙA); 


ii) La (1) è trattabile è L'unica soluzione di (6) con condizione ini- 


ziale consistente x(0) = x, è 
Da “De 
(9) x(t) = exp(-ABH)A A x + x(t). 


Anche in questo enunciato si sono seguite le notazioni prece- 
denti, ponendo di più f = (xA48)!f; ancora, k, n, Ai, APa 69, W so- 
no indipendenti da x. 

Nella (9) è evidente l'analogia con la (2): si noti che il ter- 
mine r(t) che non compare nel caso classico, scompare per A invertibile 


perché APA = TL 


Menzioniamo infine a questo punto il caso di sistemi rettan- 
i nxm diiuma TNDRTI 
golari (A, Be C con n * m), sovradeterminati o sottodeterminati: es 


Si si possono trattare secondo questa falsariga, se alla condizione (7) 


si sostituisce una ipotesi di iniettività (rispettivamente, surietti- 


vità) per una matrice (XA + B); rimandiamo a [2] o a [3]. 
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